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The origins of the Bronze Age Minoan and Mycenaean cultures 
have puzzled archaeologists for more than a century. We have 
assembled genome-wide data from 19 ancient individuals, including 
Minoans from Crete, Mycenaeans from mainland Greece, and their 
eastern neighbours from southwestern Anatolia. Here we show that 
Minoans and Mycenaeans were genetically similar, having at least 
three-quarters of their ancestry from the first Neolithic farmers 
of western Anatolia and the Aegean 1,2 , and most of the remainder 
from ancient populations related to those of the Caucasus 3 and 
Iran 4,5 . However, the Mycenaeans differed from Minoans in 
deriving additional ancestry from an ultimate source related to the 
hunter-gatherers of eastern Europe and Siberia 6-8 , introduced via 
a proximal source related to the inhabitants of either the Eurasian 
steppe 1,6,9 or Armenia 4,9 . Modern Greeks resemble the Mycenaeans, 
but with some additional dilution of the Early Neolithic ancestry. 
Our results support the idea of continuity but not isolation in the 
history of populations of the Aegean, before and after the time of 
its earliest civilizations. 

Ancient DNA research has traced the principal ancestors of early 
European farmers to highly similar Neolithic populations of Greece and 
western Anatolia, beginning in the seventh millennium bc (refs 1, 2); 
however, the later history of these regions down to the Bronze Age, 
a transformational period in the history of Eurasia 4,6,9 , is less clear. 
There is limited genetic evidence suggesting migrations from both the 
east (the area of Iran and the Caucasus), reaching Anatolia by at least 
~3800 bc (ref. 4), and the north (eastern Europe and Siberia) 
contributing Ancient North Eurasian ancestry 6,10 to all modern 
Europeans. The timing and impact of these migrations in the Aegean 
is, however, unknown. 

During the Bronze Age, two prominent archaeological cultures 
emerged in the Aegean. The culture of the island of Crete, sometimes 
referred to as ‘Minoan n , was Europe’s first literate civilization, and 
has been described as ‘Europe’s first major experience of civilization’ 12 . 
However, the Linear A syllabic ideographic and Cretan hieroglyphic 
scripts used by this culture remain undeciphered, obscuring its origins. 
Equally important was the civilization of the ‘Mycenaean culture of 


mainland Greece, whose language, written in the Linear B script, was 
an early form of Greek 13 . Cretan influence in mainland Greece and 
the later Mycenaean occupation of Crete link these two archaeological 
cultures, but the degree of genetic affinity between mainland and Cretan 
populations is unknown. Greek is related to other Indo-European 
languages, leading to diverse theories tracing its earliest speakers from 
the seventh millennium down to ^1600 bc, and proposing varying 
degrees of population change (Supplementary Information section 1). 

Genome-wide ancient DNA data provide a new source of infor¬ 
mation about the people of the Bronze Age, who were first known 
through the ancient poetic and historical traditions starting with 
Homer and Herodotus, later through the disciplines of archaeology 
and linguistics, and, more recently, by the limited information from 
ancient mitochondrial DNA 14,15 . Here we answer several questions. 
First, do the labels ‘Minoan’ and ‘Mycenaean correspond to genetically 
coherent populations or do they obscure a more complex structure of 
the peoples who inhabited Crete and mainland Greece at this time? 
Second, how were the two groups related to each other, to their neigh¬ 
bours across the Aegean in Anatolia, and to other ancient populations 
from Europe 1,2,6,8-10 and the Near East 2-5,9,16,17 ? Third, can inferences 
about their ancestral origins inform debates about the origins of their 
cultures? Fourth, how are the Minoans and Mycenaeans related to 
Modern Greeks, who inhabit the same area today? 

We generated genome-wide data from 19 ancient individuals 
(Fig. la, Extended Data Table 1 and Supplementary Information 
section 1). These comprised ten Minoans from Crete (approximately 
2900-1700 bc; labelled Minoan_Odigitria, from Moni Odigitria 
near the southern coast of central Crete; and Minoan_Lasithi, from 
the cave of Hagios Charalambos in the highland plain of Lasithi 
in east Crete). Four Mycenaeans were included from mainland 
Greece (approximately 1700-1200 bc; from the western coast 
of the Peloponnese, from Argolis, and the island of Salamis). An 
additional individual from Armenoi in western Crete (approxi¬ 
mately 1370-1340 bc; labelled Crete_Armenoi) postdated the 
appearance of Mycenaean culture on the island. Our dataset also 
included a Neolithic sample from Alepotrypa Cave at Diros Bay in 
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Figure 1 | Samples and PCA. a, Geographical locations of newly reported 
ancient data. Lines point to sampling locations; jitter is added to show 
the number of sampled individuals per location, b, Three hundred and 
thirty-four ancient individuals projected onto the first 2 principal 


components computed on a sample of 1,029 present-day West Eurasians 4,5,10,31 , 
including 30 Modern Greek samples from Greece and Cyprus. For 
abbreviations of population names, see Methods. 


the southern Peloponnese (about 5400 bc), adding to previously 
published samples from northern Greece 2 (collectively labelled 
Greece_N). Finally, it included three Bronze Age individuals (approxi¬ 
mately 2800-1800 bc; labelled Anatolia_BA) from Harmanoren 
Gondiirle in southwestern Anatolia (Turkey), adding knowledge 
about genetic variation in Anatolia after the Neolithic/Chalcolithic 
periods 1,2,4,17 (Supplementary Information section 1). We processed 
the ancient remains, extracted DNA, and prepared Illumina libraries 
in dedicated clean rooms (Methods and Supplementary Table 1), 
and, after initial screening for mitochondrial DNA, used in-solution 
hybridization 18 to capture ~ 1.2 million single nucleotide polymor¬ 
phisms (SNPs) 6,19 on the ancient samples. We assessed contamination 
by examining the rate at which they matched the mitochondrial con¬ 
sensus sequence (Supplementary Table 2) and the rate at which male 
samples were heterozygous on the X chromosome (Methods). We com¬ 
bined the dataset of the 19 ancient individuals with 332 other ancient 


individuals from the literature, 2,614 present-day humans genotyped 
on the Human Origins array, and 2 present-day Cretans (Methods). 

We performed principal component analysis (PCA) 20 (Methods), 
projecting ancient samples onto the first two principal components 
inferred from present-day West Eurasian populations 10 that form 
two south-north parallel dines in Europe and the Near East along 
principal component 2. Minoans and Mycenaeans were centrally 
positioned in the PCA (Fig. lb), framed to the left by ancient pop¬ 
ulations from mainland Europe and the Eurasian steppe, to the right 
by ancient populations from the Caucasus and Western Asia, and 
to the bottom by Early/Middle Neolithic farmers from Europe and 
Anatolia. The Neolithic samples from Greece clustered with these 
farmers and were distinct from the Minoans and Mycenaeans. The 
Bronze Age individuals from southwestern Anatolia were also distinct, 
intermediate between Anatolian and Levantine populations towards 
the bottom, and populations from Armenia, Iran, and the Caucasus 
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Table 1 | Admixture modelling of Bronze Age populations 

Ancestral sources Mixture proportions Standard errors 
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For each test population, mixture proportions from four source populations with their standard errors are given. Ancestry is inferred from both ‘ultimate’ sources representing the earliest populations, 
and ‘proximate’ sources representing populations down to the Bronze Age (Supplementary Information section 2). Column A lists ‘northern’ sources from eastern Europe and Siberia, including the 
Eurasian steppe; column B lists ‘eastern’ sources from Iran, the Caucasus, and Anatolia (after the Early Neolithic); column C lists ‘local’ sources from Anatolia and the Aegean; column D lists sources 
from the Levant. For abbreviations of population names, see Methods. 


towards the top. ADMIXTURE analysis (Methods and Extended 
Data Fig. 1) showed that Minoans and Mycenaeans both possessed 
a pink’ genetic component {K=S and greater) shared with Bronze 
Age southwestern Anatolians, Neolithic Central Anatolians from 
Tepecik-Qftlik 17 , a Chalcolithic northwestern Anatolian 1 , and western 
Anatolians from Kumtepe 16 . This component was maximized in the 
Mesolithic/Neolithic samples from Iran 4,5 and hunter-gatherers from 
the Caucasus 3 (Extended Data Fig. 1). It was not found in the Neolithic 
of northwestern Anatolia, Greece, or the Early/Middle Neolithic pop¬ 
ulations of the rest of Europe, only appearing in the populations of the 
Late Neolithic/Bronze Age in mainland Europe 6 , introduced there by 
migration from the Eurasian steppe 1,6 . 

Beyond the visual impressions of PCA and ADMIXTURE, we for¬ 
mally tested the relationships among populations from our study and 
the literature, using/ 4 -statistics of the form / 4 (X, Y ; Test, Chimp) that 
evaluated whether Test shared more alleles with X or Y. We found that 
Test populations from Iran, the Caucasus, and eastern Europe shared 
more alleles with Minoans and Mycenaeans than with the Neolithic 
population of Greece (Extended Data Fig. 2a, b). The Minoans from 
the Lasithi plateau in the highlands of eastern Crete and from the coast 
of southern Crete (Extended Data Fig. 2c) were consistent with being 
a homogeneous population. Mycenaeans differed from these Minoans 
in sharing significantly fewer alleles with Neolithic people from the 
Levant, Anatolia, Greece, and mainland Europe (Extended Data 
Fig. 2d). In comparison, the Bronze Age Anatolians shared fewer alleles 
with ancient Europeans and more with ancient populations of Iran and 
the Levant (Extended Data Fig. 3). We used/ 3 -statistics of the form 
/ 3 (Refi, Ref 2 ; Test) that, if negative, showed that Test was admixed 
from sources related to the Ref b Ref 2 source populations. We did not 
find significantly negative (Refi, Ref 2 ) pairs for Minoans or Bronze 
Age Anatolians (z> — 2.5), but did for Mycenaeans (—4.9<z<—3.0; 
Extended Data Fig. 4), involving early farmers from the Levant, 
Anatolia, Greece, and the rest of Europe as one source, and Iran or the 
Eurasian steppe or steppe-influenced Europeans as the other. 

We modelled Bronze Age populations using the qpAdm/qpWave 6 
framework (Methods and Supplementary Information section 2), 
which relates a set of ‘left’ populations (admixed population and 
ancestral source populations) with a set of right’ populations (diverse 
outgroups) and allows testing for the number of streams of ancestry 
from ‘right’ to ‘left’ and estimation of admixture proportions. This 
analysis showed that all Bronze Age populations from the Aegean and 


Anatolia are consistent with deriving most (approximately 62-86%) of 
their ancestry from an Anatolian Neolithic-related population (Table 1). 
However, they also had a component (approximately 9-32%) of 
eastern’ (Caucasus/Iran-related) ancestry. It was previously shown that 
this type of ancestry was introduced into mainland Europe via Bronze 
Age pastoralists from the Eurasian steppe, who were a mix of both 
eastern European hunter-gatherers and populations from the Caucasus 
and Iran 4,6 ; our results show that it also arrived on its own, at least 
in the Minoans, without eastern European hunter-gatherer ancestry. 
This ancestry need not have arrived from regions east of Anatolia, 
as it was already present during the Neolithic in central Anatolia at 
Tepecik-Qiftlik 17 (Supplementary Information section 2). The eastern 
influence in the Bronze Age populations from Greece and southwestern 
Anatolia is also supported by an analysis of their Y chromosomes. 
Four out of five males belonging to Minoans, Mycenaeans, and 
southwestern Anatolians (Supplementary Information section 3) 
belonged to haplogroup J, which was rare or non-existent in earlier 
populations from Greece and western Anatolia who were dominated by 
Y-chromosome haplogroup G2 (refs 1,2,17). Haplogroup J was present 
in Caucasus hunter-gatherers 3 and a Mesolithic individual from Iran 4 , 
and its spread westwards may have accompanied the eastern’ genome¬ 
wide influence. 

The Minoans could be modelled as a mixture of the Anatolia 
Neolithic-related substratum with additional ‘eastern’ ancestry, but 
the other two groups had additional ancestry: the Mycenaeans had 
approximately 4-16% ancestry from a ‘northern’ ultimate source related 
to the hunter-gatherers of eastern Europe and Siberia (Table 1), while 
the Bronze Age southwestern Anatolians may have had ~6% ancestry 
related to Neolithic Levantine populations. The elite Mycenaean 
individual from the ‘royal’ tomb at Peristeria in the western Peloponnese 
did not differ genetically from the other three Mycenaean individuals 
buried in common graves. To identify more proximate sources of 
the distinctive eastern European/north Eurasian-related ancestry 
in Mycenaeans, we included later populations as candidate sources 
(Supplementary Information section 2), and could model Mycenaeans 
as a mixture of the Anatolian Neolithic and Chalcolithic-to-Bronze Age 
populations from Armenia (Table 1). Populations from Armenia pos¬ 
sessed some ancestry related to eastern European hunter-gatherers 4 , 
so they, or similar unsampled populations of western Asia, could have 
contributed it to populations of the Aegean. This model makes geo¬ 
graphical sense, since a population movement from the vicinity of 
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Figure 2 | Genetic differentiation of Bronze Age populations to 
present-day populations. The F S t inbreeding coefficient (Methods) 
between newly reported populations and present-day West Eurasian 
populations. This shows a pattern of genetic affinity between Bronze 
Age and present-day populations from the corresponding broad 


geographical regions: a, Mycenaeans; b, Minoans from Hagios 
Charalambos (Lasithi regional unit); c, Minoans from Moni Odigitria 
(Herakleion regional unit); d, southwestern Bronze Age (BA) Anatolians. 
The same pattern also applies to Bronze Age populations from other 
regions of West Eurasia (Extended Data Fig. 5). 


Armenia could have admixed with Anatolian Neolithic-related farm¬ 
ers on either side of the Aegean. However, Mycenaeans can also be 
modelled as a mixture of Minoans and Bronze Age steppe populations 
(Table 1 and Supplementary Information section 2), suggesting that, 
alternatively, eastern ancestry arrived in both Crete and mainland 
Greece, followed by about 13-18% admixture with a northern steppe 
population in mainland Greece only Such a scenario is also plausible: 
first, it provides a genetic correlate for the distribution of shared top- 
onyms in Crete, mainland Greece, and Anatolia discovered in ref. 21; 
second, it postulates a single migration from the east; third, it proposes 
some gene flow from geographically contiguous areas to the north 
where steppe ancestry was present since at least the mid-third millen¬ 
nium bc (refs 6, 9). We validated inferences from qpAdm by treating 
source populations as ghosts’ and re-estimating mixture proportions 4 , 
by examining the correspondence between qpAdm estimates and PCA 4 
(Extended Data Fig. 5), and by comparing simulated individuals of 
known ancestry against the Mycenaeans (Extended Data Fig. 6). 

Geographical structure may have prevented the spread of the north¬ 
ern’ ancestry from the mainland to Crete, contributing to genetic 
differentiation. Such a structure may in principle, be long-standing, 
even before the advent of the Neolithic in the seventh millennium bc. 
Alternatively both northern and eastern ancestry may have arrived 
in the Aegean at any time between the Early Neolithic and the Late 
Bronze Age. Wider geographical and temporal sampling of pre- 
Bronze Age populations of the Aegean may better trace the advent 
of northern’ and eastern ancestry in the region. However, sampled 
Neolithic samples from Greece, down to the Final Neolithic ~4100 bc 
(ref. 2), do not possess either type of ancestry suggesting that the 
admixture we detect probably occurred during the fourth to second 
millennium bc time window. Other proposed migrations, such as 
settlement by Egyptian or Phoenician colonists 22 , are not discernible 


in our data, as there is no measurable Levantine or African influence 
in the Minoans and Mycenaeans, thus rejecting the hypothesis that 
the cultures of the Aegean were seeded by migrants from the old 
civilizations of these regions. On the other hand, migrants from areas 
east or north of the Aegean, while numerically less influential than the 
locals, may have contributed to the emergence of the third to second 
millennium bc Bronze Age cultures as creative disruptors’ of local tra¬ 
ditions, bearers of innovations, or through cultural interaction with 
the locals, coinciding with the genetic process of admixture 23 . Relative 
ancestral contributions do not determine the relative roles in the rise 
of civilization of the different ancestral populations; nonetheless, the 
strong persistence of the Neolithic substratum does suggest a key role 
for the locals in this process. 

Phenotype prediction from genetic data has enabled the recon¬ 
struction of the appearance of ancient Europeans 1,24 who left no visual 
record of their pigmentation. By contrast, the appearance of the Bronze 
Age people of the Aegean has been preserved in colourful frescos and 
pottery, depicting people with mostly dark hair and eyes 25 . We used the 
HIrisPlex 26 tool (Supplementary Information section 4) to infer that the 
appearance of our ancient samples matched the visual representations 
(Extended Data Table 2), suggesting that art of this period reproduced 
phenotypes naturalistically. 

We estimated the fixation index, Fst> of Bronze Age populations 
with present-day West Eurasians, finding that Mycenaeans were least 
differentiated from populations from Greece, Cyprus, Albania, and 
Italy (Fig. 2), part of a general pattern in which Bronze Age popula¬ 
tions broadly resembled present-day inhabitants from the same region 
(Extended Data Fig. 7). Modern Greeks occupy the intermediate space 
of the PCA along principal component 1 (Fig. lb) between ancient 
European and Near Eastern populations, such as those of the Bronze 
Age. They are not, however, identical to Bronze Age populations, as 
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they are above them along principal component 2 (Fig. lb). This is 
because Neolithic farmers shared fewer alleles with Modern Greeks 
than with Mycenaeans (Extended Data Fig. 8), consistent with addi¬ 
tional later admixture 27,28 . 

The Minoans and Mycenaeans, sampled from different sites in Crete 
and mainland Greece, were homogeneous, supporting the genetic 
coherency of these two groups. Differences between them were modest, 
viewed against their broad overall similarity to each other and to 
the southwestern Anatolians, sharing in both the ‘local’ Anatolian 
Neolithic-like farmer ancestry and the ‘eastern’ Caucasus-related 
admixture. Two key questions remain to be addressed by future 
studies. First, when did the common eastern’ ancestry of both Minoans 
and Mycenaeans arrive in the Aegean? Second, is the ‘northern’ ancestry 
in Mycenaeans due to sporadic infiltration of Greece, or to a rapid 
migration as in Central Europe 6 ? Such a migration would support the 
idea that proto-Greek speakers 29 formed the southern wing of a steppe 
intrusion of Indo-European speakers. Yet, the absence of‘northern 
ancestry in the Bronze Age samples from Pisidia, where Indo-European 
languages were attested in antiquity, casts doubt on this genetic- 
linguistic association, with further sampling of ancient Anatolian 
speakers needed. Whatever the answer to these questions, the discovery 
of at least two migration events into the Aegean in addition to the first 
farming dispersal before the Bronze Age, and of additional popula¬ 
tion change since that time, supports the view that the Greeks did not 
emerge fully formed from the depths of prehistory, but were, indeed, a 
people ever in the process of becoming’ 30 . 

Online Content Methods, along with any additional Extended Data display items and 
Source Data, are available in the online version of the paper; references unique to 
these sections appear only in the online paper. 
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METHODS 

No statistical methods were used to predetermine sample size. The experiments 
were not randomized and the investigators were not blinded to allocation during 
experiments and outcome assessment. 

Ancient DNA. An overview of which steps in processing the ancient samples were 
undertaken in which laboratory is provided in Supplementary Table 1. 

Dublin, Ireland. The inner ear area of each petrous bone was identified, isolated, 
then ground to a fine powder. Cleaning and isolation of the cochlea was performed 
using aluminium oxide powder in a sandblasting chamber. Once isolated, it was 
decontaminated by ultraviolet irradiation for 7.5 min on each side, ground on a 
mixer mill to a weight of about 50 mg, and finally transferred to a sterile Eppendorf 
tube. All procedures were conducted in clean and dedicated ancient DNA facilities. 
Seattle, Washington, USA. Teeth processed in this laboratory were decontam¬ 
inated and pulverized to powder in clean and dedicated ancient DNA facilities 
following previously published methods 11 . 

Leipzig, Germany. As previously described 32 , sampling, extraction, and prepa¬ 
ration of double-indexed, double-stranded libraries took place in the clean room 
facilities of the Max Planck Institute for Evolutionary Anthropology, Leipzig, 
Germany (MPI-EVA), followed by enrichment of human mitochondrial DNA 33 . 
Enriched libraries were sequenced on an Illumina GAIIx platform for 2 x 76 + 7 
cycles and the resulting data were mapped to the revised Cambridge Reference 
Sequence using the EAGER pipeline to evaluate DNA preservation (Supplementary 
Table 2). These libraries were then shipped to Boston, Massachusetts, USA, where 
nuclear target enrichment was performed (see below). 

Tubingen, Germany. Pre-PCR steps took place in the clean room facilities of 
the Institute for Archaeological Sciences at the University of Tubingen, Germany. 
After surface irradiation with ultraviolet light, the tooth was sawn apart trans- 
versally at the border of crown and root, and dentine powder from the inside 
the crown was sampled using a sterile dentistry drill. Extraction, library prepara¬ 
tion, and enrichment of human mitochondrial DNA used the same protocols as 
described for MPI-EVA, with the addition of an updated extraction protocol 34 . 
Sequencing of shotgun and mitochondrial-DNA-enriched libraries took place at 
the facilities of the Frauenklinik of the University of Tubingen, on an Illumina 
MiSeq for 2 x 150 + 8 cycles or on an Illumina HiSeq 2500 for 2 x 101 + 8 cycles 
(Supplementary Table 2). 

Additional libraries were produced including full or partial 35 repair with 
uracil-DNA glycosylase and endonuclease VIII to remove deaminated bases. 
In-solution enrichment was performed using previously reported protocols 6,18 . 
Two SNP sets of 394,577 SNPs ( c 390k capture’ 6 ) or 1,237,207 SNPs (‘1240k cap¬ 
ture’ 1 ) were targeted. Sequencing took place in the facilities of the Frauenklinik, 
University of Tubingen, on an Illumina HiSeq 2500 for 2 x 101 + 8 cycles and at 
the facilities of the University of Kiel on a HiSeq 4000 for 2 x 150 + 8 cycles. One 
uracil-DNA glycosylase-treated library (10071) was sent to Boston, Massachusetts, 
for nuclear target enrichment (see below). 

Boston, Massachusetts, USA. The bone powders, prepared from petrous bones in 
Dublin, Ireland, were sent to Boston, where DNA extractions and barcoded library 
preparations without uracil removal were performed in the Harvard Medical 
School cleanroom following previously described protocols 34-36 . At the screening 
stage, libraries were (1) shotgun sequenced and (2) sequenced after enriching for 
the human mitochondrial DNA 37 together with some nuclear loci to approximate 
the nuclear coverage and mitochondrial contamination. 

All four libraries (barcoded) prepared in Boston, three libraries (indexed) pre¬ 
pared in Leipzig, and one library (indexed) prepared in Tubingen, were used to 
perform 390k (ref. 6) and 840k (ref. 19) or 1240k (=390k + 840k) targeted capture 
of a total of 1,233,013 SNPs, following the in-solution target enrichment protocol 
in ref. 18 and sequenced either on an Illumina HiSeq 2500 or an Illumina NextSeq 
500 (see Supplementary Table 1 for details). 

For each sample, each SNP position was represented by a randomly chosen 
sequence, restricting to those with a minimum mapping quality (MAPQ > 10), 
sites with a minimum sequence quality (>20), and removing two bases at the 
ends of reads 4 . 

Testing for contamination. Modern human contamination of the mitochondrial 
DNA was assessed using the software schmutzi 38 , which took into account that 
the consensus sequence should be reconstructed from reads showing character¬ 
istics of ancient DNA and originating from a single individual (Supplementary 
Table 2). We assessed contamination by examining heterozygosity on the 
X chromosome in five males (possessing only one copy of the X chromosome) using 
ANGSD 39 (Supplementary Information section 3); this was in the range 0.3-4%. 
Indirect evidence that the females in our dataset (for which X-chromosome-based 
contamination estimation was impossible) were authentic was furnished by their 
clustering with male samples and distinctiveness from present-day Greek or central 
European populations that may have possibly contaminated them (Fig. lb). We also 


computed/ 4 -statistics of the form/ 4 (Males, Females; Test, Chimp) for populations 
that had both male and female individuals for all ancient or present-day Test pop¬ 
ulations in our dataset. If female samples were substantially contaminated from 
a source related to Test, these statistics would be significantly negative; however, 
we found that the Z-score of these statistics was —1.6 < Z< 2.5. We thus included 
both male and female samples in our analysis to maximize sample size instead of 
restricting it to damaged molecules for females 8 . 

Modern human data. We used a dataset of 2,614 individuals genotyped on the 
Affymetrix Human Origins array 4,5,10,31 , including 28 Modern Greek (from Greece 
and Cyprus) samples previously described 10 . We also included data from two 
Modern Greeks from Crete whose whole-genome sequences were published as 
part of the Simons Genome Diversity Project 40 . We also analysed Modern Greek 
data from Thessaly and Central Greece 41 and diverse regions 27,42 genotyped on 
Illumina arrays. 

Datasets. We analysed two datasets: HO, which includes the Affymetrix Human 
Origins genotyping data together with 351 ancient humans (including samples from 
the literature 1-5,7-10,16,17,43-51 and the newly reported data) on 591,642 autosomal 
SNPs; and the HOIll dataset, which does not include the Human Origins data, but 
has a larger number of 1,054,671 autosomal SNPs 4 . We did not use previously per¬ 
formed genotype calls of data from the literature, but re-processed them, beginning 
with the original data release format (FASTQ or BAM). The main analysis dataset 
was HOIll, except for analyses that included modern populations, in which case 
the HO dataset was analysed. For the analysis of Illumina genotype data of Modern 
Greeks (Extended Data Fig. 6), a total of489,148 autosomal SNPs were analysed. 
Abbreviations used. For brevity, we used the following abbreviations in population 
names, following the convention of ref. 4: CHG, Caucasus hunter-gatherers; EHG, 
Eastern European hunter-gatherers; WHG, Western European hunter-gatherers; 
SHG, Scandinavian hunter-gatherers; N, Neolithic; EN, Early Neolithic; 
MN, Middle Neolithic; ChL, Chalcolithic; LNBA, Late Neolithic/Bronze Age; 
BA, Bronze Age; EBA, Early Bronze Age; EMBA, Early/Middle Bronze Age; 
MLBA, Middle/Late Bronze Age; IA, Iron Age. 

PCA. PCA was performed in the smartpca program of EIGENSOFT 20 , using 
default parameters and the lsqproject: YES 10 and numoutlieriter: 0 options. PCA 
was performed on 1,029 present-day West Eurasians and 334 ancient samples were 
projected (Fig. lb); Upper Palaeolithic individuals before the appearance of the 
Villabruna cluster 8 plot in the middle of present-day West Eurasian variation and 
are not shown. 

ADMIXTURE analysis. ADMIXTURE analysis 52 of the HO dataset was per¬ 
formed after pruning for linkage disequilibrium in PLINK 53,54 with parameters 
indep-pairwise 200 25 0.4, after which 299,971 SNPs were retained. Twenty repli¬ 
cates of the analysis were performed with different random seeds, and the highest 
likelihood replicate for each value of K was retained. We show the K=2to K= 17 
results for the 351 ancient and 30 Modern Greek samples in Extended Data Fig. 1. 
/-statistics. The/- and/ 4 -statistics were computed in ADMIXTOOLS 31 using 
the programs qp3Pop and qpF4ratio with default parameters, and qpDstat with 
f4mode: YES. Standard errors were computed with a block jack-knife 55 . When an 
ancient population was the target for/-statistics, we set the inbreed: YES parame¬ 
ter, as our data were represented by pseudo-haploid genotypes, which introduced 
artificial genetic drift that masked the negative signal of admixture 31 . 

Testing for the number of streams of ancestry and estimating mixture propor¬ 
tions. We used the qpWave 6,56,57 /qpAdm 6 framework, which relates a set of‘left’ 
populations (the population of interest and candidate ancestral sources) to a set 
of‘right’ populations (diverse outgroups), testing for the number of streams of 
ancestry from ‘right’ to ‘left’ and estimating mixture proportions. 

Simulations of admixed individuals. We simulated admixed individuals 
(Supplementary Information section 2) given a set of sources and mixture propor¬ 
tions by first sampling (at each SNP) one of the sources (according to the mixture 
proportions), and then one of the individuals from that population (with equal 
probability). Because of missingness, the data-generating mixture proportions 
did not correspond precisely to the actual ancestry of simulated individuals and 
we corrected for this bias (Supplementary Information section 2). We noted the 
maximum absolute value of the Z-score of the statistic/(Mycenaean, Simulated; 
A, B), where A, B were two outgroup populations to test whether, for a particular 
choice of ancestry of Simulated, it formed a clade with the sampled Mycenaeans. 
Estimation of F S t coefficients. We estimated +st in smartpca 20 with the default 
parameters inbreed: YES 57 , and fstonly: YES. 

Phenotypic inference. The ancient samples had low coverage (median 0.87 x) 
and thus diploid genotypes could not be reliably assessed for them. However, we 
could use the low coverage data to compute allele frequencies in all individuals and 
the Bronze Age Aegean using a likelihood approach 1 . We then sampled from the 
posterior distribution of the genotypes g given the read counts r of the reference 
allele and t of the total reads covering a site. We took 100 random genotype samples 
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per individuals and submitted them to HIrisPlex 26 , obtaining an estimate of the 
uncertainty of phenotype inference (Supplementary Information section 4 and 
Extended Data Table 2). 

Data availability. The aligned sequences are available through the European 
Nucleotide Archive under accession number PRJEB20914. Genotype datasets used 
in analysis are available at https://reich.hms.harvard.edu/datasets. All other data 
are available from the corresponding authors upon reasonable request. 
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Extended Data Figure 1 | ADMIXTURE analysis. ADMIXTURE analysis (Methods) with K — 2 to K= 17 is shown. Three hundred and fifty-one 
ancient and 2,616 present-day individuals were used in this analysis; ancient samples and present-day Greeks are displayed. To avoid visual clutter of 
labels, individuals in populations with sample size < 5 are shown with thicker lines. 
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Extended Data Figure 2 | Symmetry testing of Aegean Bronze Age 
populations. The statistic f^X, Y; Test, Chimp) is shown with ±3 standard 
errors. Each panel is titled with the pair X, Y. Populations are ordered 
according to the value of the statistic. Positive values indicate that Test 
shares more alleles with X than Y , and negative values that it shares more 
with Y than X. a, ‘Northern and ‘eastern populations share more alleles 


-o.5io -0.605 o.doo ados 

with Minoans than with Neolithic Greece, b, ‘Northern and ‘eastern 
populations share more alleles with Mycenaeans than with Neolithic 
Greece, c, Minoans from Lasithi and Moni Odigitria are symmetrically 
related to diverse populations, d, Neolithic populations from Anatolia, 
Europe, Greece, and the Levant share fewer alleles with Mycenaeans than 
with Minoans. 
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Extended Data Figure 3 | Symmetry testing of Anatolian Bronze Age 
populations. The statistic/ 4 (X, Y ; Test, Chimp) is shown with ±3 standard 
errors. Each panel is titled with the pair X, Y. Populations are ordered 
according to the value of the statistic. Positive values indicate that Test 
shares more alleles with X than Y, and negative values that it shares more 
with Y than X. a, European, Siberian, and Caucasus hunter-gatherers share 
fewer alleles with Bronze Age Anatolians from Harmanoren Gondiirle 


-0.605 o.doo o.do5 o!3io 

than with a Chalcolithic Anatolian from Barcin. b, Bronze Age Anatolians 
differ from Neolithic ones in sharing more alleles with populations of Iran, 
the Caucasus, and the Steppe than with those of Europe, c, Bronze Age 
Anatolians differ from Minoans in sharing more alleles with populations 
from Neolithic Iran than Neolithic Anatolia and Europe, d, Bronze Age 
Anatolians differ from Mycenaeans in sharing more alleles with Neolithic 
and Bronze Age populations of the Levant. 
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Extended Data Figure 4 | The/ 3 -statistics of Mycenaeans as a target with 
different pairs of reference populations. The value of the statistic 
/ 3 (Refi, Ref 2 ; Mycenaean) with ±3 standard errors; only the population 


- 1 - 1 - 1 

- 0.01 0.00 0.01 

pairs (Refi, Ref 2 ) for which the Z-score of the statistic is less than —2 
are shown. Negative values indicate that the Mycenaean population is 
admixed from sources related to the two reference populations. 
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Extended Data Figure 5 | Correspondence of qpAdm estimates with 
PCA. As a way of validating qpAdm models of admixture for Mycenaeans 
from three ancestral populations (Anatolia_N or Minoan_Lasithi), 
(Armenia_ChL or Armenia_MLBA), (Steppe_EMBA, Steppe_MLBA, 


Europe_LNBA), representing substratum, eastern, and northern ancestry, 
respectively (Supplementary Information section 2), we plot the qpAdm- 
predicted position in the PCA space of Fig. 1 versus the actual position of 
the Mycenaean population. 


© 2017 Macmillan Publishers Limited, part of Springer Nature. All rights reserved. 






LETTER 


RESEARCH 


54%Anatolia_N,46%Armenia_ChL,0%Steppe_EMBA 
61 %Anatolia_N ,36%Armenia_ChL,3%Steppe_EMBA 
(Max difference: 10%) 


Steppe EMBA Best estimate 



o 0 o°o°o%; 

m PvO p-op;a-^-^ 

80 oOOOOOC 

O O 0 -0-0 o o 
oooooooc 
OOOooooo 
DDOOoo o O c 

OOooooo o o 

Anatolia_N 4 


qpAdm (p=0.57) 
simulation (|Z|=2.6) 


54%Anatolia_N,46%Armenia_ChL,0%Steppe_MLBA 
59%Anatolia_N,38%Armenia_ChL,4%Steppe_MLBA 
(Max difference:8%) 


Max |Z| 


O 12 
021 


Steppe_MLBA 


Best estimate 


qpAdm (p=0.59) 
simulation (|Z|=2.6) 



* 


Armenia_ChL 


0°p"QO OOOOOQC^- 

oo o o cKKkSS^aas 

O 0:0 0; OrO-O-O-OC 
80 OOOOOOOO O O 
OOOOOO OOOOC 

00O ooooooOO* 
lOO o o o o o«o o o O 
OOOOo o o o o O o < 

OOOO O o o o o-*- o-o 
Anatolia_N 40 60 Armenia_ChL 


54%Anatolia_N,46%Armenia_ChL,0%Europe_LNBA 
58%Anatolia_N,39%Armenia_ChL,4%Europe_LNBA 
(Max difference:8%) 


Max |Z| 

• 2.6 
O 5 


Europe_LNBA 

O 


Best estimate 


o o 

cggtL, 

° 7 ° 

018 40 o o 0 0_ O ^60 

018 , 0 ?g 0 o°o°o°o°oc> 

60 O 0 O%%^O°O°Q 0 9_ r 
OOO OOCbOO U - 1 
o O OlO OcO. 000 c 

o o v o o o oooof 

OOOO OOOO O Q 
000OOO 000C 

O OOOOOO 0 _Q. 

) o O o O O OOO r 
\ooo A ooo r 
°o°o° o°o OOO 

3 0^0 0 o • 0-6 

40 


♦ qpAdm (p=0.59) 

• simulation (|Z|=2.6) 



ArmeniaChL 


62%Anatolia_N,31%Armenia_MLBA,6%Steppe_EMBA 
66%Anatolia_N,31%Armenia_MLBA,3%Steppe_EMBA 
(Max difference:4%) 



62%Anatolia_N,31%Armenia_MLBA,7%Steppe_MLBA 
65%Anatolia_N,33%Armenia_MLBA,3%Steppe_MLBA 
(Max difference:4%) 



64%Anatolia_N,32%Armenia_MLBA,5%Europe_LNBA 
64%Anatolia_N,33%Armenia_MLBA,2%Europe_LNBA 
(Max difference:2%) 



78%Minoan_Lasithi,0%Armenia_ChL,22%Steppe_EMBA 
90%Minoan_Lasithi,-2%Armenia_ChL, 12%Steppe_EMBA 
(Max difference: 13%) 


Max |Z| 


Steppe.EMBA B est estimate 

• 3.3 rQQs ♦ qpAdm (P =0 - 25 ) 

OOO • simulation (|Z|=3.3) 

0 6 20 OOOOrq 

010 20 ooooc ? 0 

OOOOOO 

015 OOOOOOO 

029 40 OOOOOOOO60 

„ 000000000 

029.6 OOOOOOO OOO 

o-o 00 0000000 

60 00000000000040 
o o O OO OOOOOOOO 
O o o o o o OOOOOOOO 
^ OoooooOOOOOOOOO 

♦ I oVcfoooooooOOOOOOon 

O O O oooooooooOOOOCT 
-- 


5 o o o o o o OOOOOOOOO 
>a o-ojo'o"070707Q7070:0;0 _ 010 

o o o o o o o o o 000OOOOOOOO 
o-o-o-o-o-oo-oo-o-oo-o-o-O-O-O-O-OOO 

Minoan_Lasithi 40 60 Armenia_ChL 


61 %Minoan_Lasithi,21 %Armenia_ChL, 18%Steppe_MLBA 
81%Minoan_Lasithi,6%Armenia_ChL,12%Steppe_MLBA 
(Max difference:21%) 

Max |Z| Bas, esUmate 

•3 ♦ qpAdm (p=0.41) 

o 5 nnnn * Simulation (|Z|=3) 

09 20 OO(Wf 0 

lu o°o 0 o 0 o 0 o 0 o°o 

021 »o 0 o 0 o 0 o°om?f 

021.4 OOOOOOOOOO 

OOOOOOOOOOO 
60 OOOOOOOOOOOO 40 
O o 00 OOOOOOOOO 
000000 O OOOOOOO 
o o o o o 00 OOOOOOOO 
<fo o o o oooooooOOOO ? n 

80 oooo*ooooOO OOOOOCr 
ooooooooooOOOOOOOO 
00000000000 OOOOOOOO 
00000000000 OOOOOOOOO 
o o 0 0 0 o 0 0 6-0 0 0 O O OOOOOOO 
Minoan_Lasithi 40 60 Armenia_ChL 


79%Minoan_Lasithi,0%Armenia_ChL,21%Europe_LNBA 
78%Minoan_Lasithi,10%Armenia_ChL,12%Europe_LNBA 
(Max difference: 10%) 


Max |Z| 


Europe_LNBA 

O 


Best estimate 


• 2.7 O O 

04 oXxb .. 

07 20 oood 

o 10 OOO 0(0 

00 mmb 
018 40 o o o o o o (f 3 
018.1 OOOOOOOO 

0000.000 010 

60 OOO OOO O OOO40 
bU OOOOOOOO'"' 
OOOOOOOO 
o 0000OOOC 
0*0 00000OO 
80 • o o o o o o O O C 
000000000O 
0000000OOO'" 


♦ qpAdm (p=0.41) 

• simulation (|Z|=2.7) 


_ - - o-'o o 0 0 OOO 
OOO 000000 OOC. 

Minoan_Lasithi 40 60 



80%Minoan_Lasithi,2%Armenia_MLBA,18%Steppe_EMBA 
88%Minoan_Lasithi,1%Armenia_MLBA,11%Steppe_EMBA 
(Max difference:9%) 



71%Minoan_Lasithi,8%Armenia_MLBA,21%Steppe_MLBA 
82%Minoan_Lasithi,6%Armenia_MLBA, 12%Steppe_MLBA 
(Max difference: 12%) 



72%Minoan_Lasithi,12%Armenia_MLBA,16%Europe_LNBA 
80%Minoan_Lasithi,8%Armenia_MLBA,12%Europe_LNBA 
(Max difference:7%) 



Extended Data Figure 6 | Comparison of Mycenaeans and simulated 
admixed populations. We simulate admixed individuals with known 
ancestry from three ancestral populations (Anatolia_N or Minoan_ 
Lasithi), (Armenia_ChL or Armenia_MLBA), (Steppe_EMBA, Steppe_ 
MLB A, Europe_LNBA), representing substratum, eastern, and ‘northern 
ancestry, respectively (Methods and Supplementary Information section 2). 


The maximum \Z \-score of statistics/ 4 (Mycenaean, Simulated; Outgroupi, 
Outgroup 2 ) is plotted with circles of varying size (proportional to log(|Z|)) 
for each assignment of ancestry proportions. The best estimate (red) 
corresponds to the proportions that minimize \Z\, and they are compared 
against the qpAdm estimate for the same ancestral sources (blue). 
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Extended Data Figure 7 | F s t between Bronze Age and present-day 
West Eurasian populations, a, The population of Early Bronze Age 
Armenia 4 shows an affinity to present-day populations from Armenia, 
Anatolia, the Caucasus, and Iran, as does (b) Middle/Late Bronze Age 
Armenia 4,9 , c, The Bronze Age Levant 4 has an affinity to Levantine and 
Arabian populations, d, Late Neolithic/Bronze Age Europeans 1,6,9,43 most 
resemble present-day northern/central Europeans, as do (e) Early/Middle 


0 Steppe_MLBA 





y 

0.008 


Bronze Age steppe populations 1,6,9 , who also resemble populations of the 
northeast Caucasus, while (f) Middle/Late Bronze Age steppe populations 
resemble central/northern Europeans 1,9 . Jewish populations are plotted 
with a square to distinguish them from non-Jewish populations from the 
same geographical area. The plots for the newly reported populations of 
Mycenaeans, Minoans, and Bronze Age Anatolians are shown in Fig. 2. 
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Extended Data Figure 8 | Symmetry testing of Mycenaeans with Modern 
Greek populations. The statistic/^Mycenaean, Modern Greek; Test, 
Chimp) is shown with ±3 standard errors. Modern Greeks share fewer 
alleles with Levantine/Anatolian/European Neolithic populations and 
with Minoans than Mycenaeans, suggesting a dilution of Early Neolithic 


ancestry since the Bronze Age. Human Origins genotype data: a, Greeks 
from the Coriell repository 10 ; b, Greeks from Thessaloniki 10 ; c, Cypriots 10 . 
Whole-genome data: d, Cretans 40 . Illumina genotype data: e, Greeks from 
Thessaly 41 ; f, Greeks from central Greece 41 ; g, Greeks from the study in 
ref. 27. 
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Extended Data Table 1 | Information on ancient samples reported in this study 


IndividualJD 

GenotypeJD 

OtherJD 

Source 

Date 

Population_Label 

Location 

Country 

Latitude 

Longitude 

Sex 

Coverage 

Autosomal_SNPs 

mtDNA 

Y-chromosome 

12937 

12937 

A2197 

1240K 

5419±41 cal BC 

Greece_N 

Diros, Alepotrypa Cave 

Greece 

36.64 

22.38 

F 

0.870 

481848 

K1a26 


10071 

10071 

Lasithi4 

1240K 

2000-1700 BCE 

Minoan Lasithi 

Hagios Charalambos Cave, Lasithi, Crete 

Greece 

35.08 

25.83 

F 

7.312 

953157 

U5a1 


10070 

10070 

Lasithi2 

1240K 

2000-1700 BCE 

Minoan Lasithi 

Hagios Charalambos Cave, Lasithi, Crete 

Greece 

35.08 

25.83 

M 

1.267 

619767 

H13a1 

J2a1d 

10073 

10073 

Lasithi7 

1240K 

2000-1700 BCE 

Minoan Lasithi 

Hagios Charalambos Cave, Lasithi, Crete 

Greece 

35.08 

25.83 

M 

1.481 

643360 

H 

J2a1 

10074 

10074 

Lasithi9 

1240K 

2000-1700 BCE 

Minoan Lasithi 

Hagios Charalambos Cave, Lasithi, Crete 

Greece 

35.08 

25.83 

F 

0.874 

506434 

H5 


19005 

19005 

Lasithil 7 

1240K 

2000-1700 BCE 

Minoan Lasithi 

Hagios Charalambos Cave, Lasithi, Crete 

Greece 

35.08 

25.83 

F 

1.351 

388859 

H 


19006 

19006 

Salamis31 

1240K 

1411-1262 cal BCE 

(3067 ± 25 BP, DEM-2905) 

Mycenaean 

Agia Kyriaki, Salamis 

Greece 

37.97 

23.50 

F 

1.387 

361193 

X2d 


19123 

19123 

S-EVA 1263 

Armenoi 503 

1240K 

1370-1340 BCE 

Crete_Armenoi 

Armenoi, Crete 

Greece 

35.45 

24.17 

F 

0.041 

45158 

U5a1 


19127 

19127 

12V t2 

1240K 

2900-1900 BCE 

Minoan Odigitria 

Moni Odigitria, Heraklion, Crete 

Greece 

35.05 

24.81 

F 

0.035 

36475 

J2b1a1 


19128 

19128 

13V t2 

1240K 

2900-1900 BCE 

Minoan Odigitria 

Moni Odigitria, Heraklion, Crete 

Greece 

35.05 

24.81 

F 

0.016 

17081 

15 


19129 

19129 

14V t2 

1240K 

2900-1900 BCE 

Minoan Odigitria 

Moni Odigitria, Heraklion, Crete 

Greece 

35.05 

24.81 

F 

0.063 

63986 

H+163 


19130 

19130 

16V Tholos 

1240K 

2900-1900 BCE 

Minoan Odigitria 

Moni Odigitria, Heraklion, Crete 

Greece 

35.05 

24.81 

M 

0.086 

92186 

U3b3 

G2a2b2 

19131 

19131 

19V t2 

1240K 

2900-1900 BCE 

Minoan Odigitria 

Moni Odigitria, Heraklion, Crete 

Greece 

35.05 

24.81 

F 

0.095 

96946 

K1a2 


19010 

19010 

Galatas19 

1240K 

1700-1200 BCE 

Mycenaean 

Galatas Apatheia, Peloponnese 

Greece 

37.50 

23.45 

F 

0.379 

242265 

X2 


19033 

19033 

Peristeria4 

1240K 

1416-1280 cal BCE 
(3084 ± 24 BP, DEM-2903) 

Mycenaean 

Peristeria Tryfilia, Peloponnese 

Greece 

36.92 

21.70 

F 

0.439 

248912 

H 


19041 

19041 

Galatas4 

1240K 

1700-1200 BCE 

Mycenaean 

Galatas Apatheia, Peloponnese 

Greece 

37.50 

23.45 

M 

1.558 

417898 

X2 

J2a1 

12495 

12495 

A4-1 

1240K 

2558-2295 calBCE 

(3925±35 BP, Poz-81111) 

Anatolia_BA 

Harmanoren-Gondurle Hoytik, Isparta 

Turkey 

37.92 

30.71 

M 

1.981 

637146 

H 

Jla 

12499 

12499 

UC1 

1240K 

2836-2472 calBCE 
(4040±35 BP, Poz-82213) 

Anatolia_BA 

Harmanoren-Gondurle Hoytik, Isparta 

Turkey 

37.92 

30.71 

F 

0.285 

243348 

K1a2 


12683 

12683 

G3-95 

1240K 

2500-1800 BCE 

Anatolia_BA 

Harmanoren-Gondurle Hoytik, Isparta 

Turkey 

37.92 

30.71 

F 

3.695 

749308 

T2b 



Dates marked simply as bce (Before Common Era) are based on the associated archaeology of the samples. Dates marked as calBCE are based on radiocarbon dating of the samples (Supplementary 
Information section 1). 
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Extended Data Table 2 | Phenotypic inference of ancient individuals 


ID 

Q) 

>. 

Ill 

0) 

3 

00 

Q. 

Population 

PlntermediateEye 

PBrownEye 

PBIondHair 

PBrownHair 

PRedHair 

PBIackHair 

PLightHair 


PDarkHair 

Hair Color 

Eye Color 

12495 

Anatolia_BA 1.6 (4.4) 

3.6 (3.9) 

94.9 (8.3) 

10.7(6.1) 

51.6 (6.4) 

0.1 (0.1) 

37.6 (9.3) 

18.0 (11.7) 

82.0 

(11.7) 

Brown 

Brown 

12499 

Anatolia_BA 16.6 (28.3) 

7.4 (2.2) 

76.0 (287) 

2.2 (2.2) 

64.7(11.8) 

2.0 (5.3) 

31.1 (13.8) 

12.9 (20.1) 

87.1 

(20.1) 

Brown 

Blue or Brown 

12683 

Anatolia_BA 0.3 (0.9) 

1.3(1.7) 

98.4 (2.6) 

3.3 (2.5) 

33.0 (4.6) 

0.0 (0.0) 

63.7 (7.0) 

4.9 (4.5) 

95.1 

(4.5) 

Black 

Brown 

12937 

Greece_N 0.3 (1.3) 

2.2 (1.9) 

97.5 (3.2) 

3.6 (1.9) 

33.9 (6.2) 

0.1 (0.0) 

62.4 (7.4) 

6.7 (4.3) 

93.3 

(4.3) 

Black 

Brown 

10070 

Minoan_Lasithi 0.4 (1.8) 

2.2 (1.9) 

97.4 (3.7) 

30.4 (5.1) 

66.4 (5.9) 

3.2 (0.9) 

0.0 (0.0) 

100.0 (0.0) 

0.0 

(0.0) 

Brown 

Brown 

10071 

Minoan_Lasithi 0.0 (0.0) 

0.2 (0.0) 

99.8 (0.0) 

0.4 (0.0) 

20.3 (0.0) 

0.0 (0.0) 

79.3 (0.0) 

0.5 (0.0) 

99.5 

(0.0) 

Black 

Brown 

10073 

Minoan_Lasithi 0.1 (0.7) 

17(1.4) 

98.2 (2.2) 

12.5(3.4) 

61.1 (1.2) 

0.2 (0.1) 

26.2 (2.7) 

32.4 (8.8) 

67.6 

(8.8) 

Brown 

Brown 

10074 

Minoan_Lasithi 0.0 (0.0) 

1.3 (0.3) 

98.7 (0.4) 

9.3 (3.2) 

54.8 (8.5) 

0.1 (0.1) 

35.8 (10.5) 

18.8 (10.3) 

81.2 

(10.3) 

Brown 

Brown 

19005 

Minoan_Lasithi 5.2 (0.0) 

11.6 (0.0) 

83.2 (0.0) 

49.6(1.4) 

38.8(1.2) 

4.2 (0.5) 

7.4 (0.7) 

85.6 (1.7) 

14.4 

(1.7) 

Blond or Brown 

Brown 

19006 

Mycenaean 0.0 (0.0) 

1.1 (0.4) 

98.9 (0.4) 

8.7 (4.9) 

59.9 (6.4) 

1.8 (2.9) 

29.6 (11.8) 

25.7 (16.5) 

74.3 

(16.5) 

Brown 

Brown 

19033 

Mycenaean 0.4 (1.0) 

1.6 (1.9) 

98.0 (3.0) 

4.6 (3.9) 

51.0 (6.3) 

0.1 (0.5) 

44.2 (9.8) 

10.5 (13.2) 

89.5 

(13.2) 

Brown 

Brown 

19041 

Mycenaean 1.4 (0.5) 

5.3 (1.0) 

93.3(1.4) 

7.8 (0.7) 

63.2 (2.0) 

0.2 (0.4) 

28.7 (2.3) 

21.2 (2.5) 

78.8 

(2.5) 

Brown 

Brown 


We list the probability assignments for different phenotypes by HlrisPlex 26 and an assessment of the phenotype. We generate 100 random replicates of the genotypes of each individual, listing the 
standard deviation in parentheses (Supplementary Information section 4). 
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